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We study the class of two-field inflationary universe models A^ 4 /4 + g 2 X 2( t )2 /2, in which paramet- 
ric resonance during the initial stages of reheating can lead to an exponential amplification of the 
amplitude of cosmological fluctuations. Employing both analytical arguments and numerical simu- 
lations, we determine the time at which backreaction of fluctuations on the background fields shuts 
off the exponential growth, making use of the Hartree approximation, and including scalar metric 
perturbations. For the case g 2 /A = 2, we find that the amplitude of fluctuations after preheating will 
exceed the observational upper bound independent of the value of A, unless the duration of inflation 
is very long. Cosmological fluctuations are acceptably small for g 2 /X > 8. We also find that the 
addition of x-field self-interaction can limit the growth of fluctuations, and in the negative-coupling 
case the system can become effectively single-field, removing the resonance. 
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I. INTRODUCTION 

It has long been realized that reheating is a crucial part 
of the inflationary scenario. During reheating the large 
energy density contained within the coherently oscillat- 
ing inflaton field is converted into particle excitations of 
whatever fields are coupled to the inflaton, vastly increas- 
ing the temperature and entropy density and setting the 
stage for the standard big bang phase. If inflation is ever 
to be a useful picture for describing the early universe, 
then it is essential to understand the details of how the 
vacuum energy is transformed into familiar particles. 

Reheating can occur very efficiently through the pro- 
cess of parametric resonance |l]-[| . Field modes within 
certain resonance bands in fc-space grow exponentially 
with time, defining the "preheating" era. The possibility 
of resonant growth of linear scalar metric perturbations 
was first studied in Q . Recently it has been argued that 
the resonance of scalar metric perturbations can extend 
to k <C aH, i.e. that super-Hubble perturbations can be 
amplified J^|| . This opens up the possibility of new ob- 
servational consequences, since the scales relevant to the 
cosmic microwave background and large-scale structure 
are much larger than the Hubble radius during preheat- 
ing. The importance of the gauge-invariant formalism 
for cosmological perturbations and the study of the 
"traditionally conserved" Bardeen parameter £ was em- 
phasized in where it was found that simple single- 
field chaotic inflation models do not exhibit super-Hubble 
growth beyond what is expected in the absence of para- 
metric resonance. The absence of parametric amplifica- 
tion of super-Hubble modes in these single field models 
was shown to hold in a full nonlinear treatment ||, and 
a general no-go theorem in these models was suggested 
in @. 

For the first model which was claimed to exhibit 



growth of super-Hubble metric perturbations beyond 
that of the usual theory of reheating ^ (see also 
it was soon realized that the growth was unimportant 
since it followed a period of exponential damping during 
inflation JEi|-p"4|| . This damping of super-Hubble modes 
arises because the field perturbations which are amplified 
during preheating have an effective mass greater than the 
Hubble parameter during inflation. This results in a very 
"blue" power spectrum at the end of inflation, with a se- 
vere deficit at the largest scales [[TiJ . The relatively plen- 
tiful small-scale modes can also grow resonantly during 
preheating. Thus the end of parametric resonance oc- 
curs when the backreaction of the dominant small-scale 
modes becomes important, and the cosmological-scale 
modes are still negligible. An obvious class of models 
to study, then, consists of those with small masses dur- 
ing inflation and strong super-Hubble resonance |^5| , |l6| . 
A simple example was provided by Bassett and Viniegra 
||17|| , namely that of a massless self-coupled inflaton <p 
coupled to another scalar field Xi i- e - a model with po- 
tential 
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This model has been studied in detail, but in the absence 
of metric perturbations, by Greene et al. fl8|| , who found 
that the model contains a strong resonance band for \ 
fluctuations which extends to k = for the choice g 2 — 
2A. Bassett and Viniegra (l7| found that super-Hubble 
metric perturbations arc resonantly amplified as well in 
this model (see also p6[). 

To date, however, none of the linearized analyses of 
parametric amplification of super-Hubble-scale metric 
fluctuations in the model (Jf]) has included the effects of 
backreaction on the evolution of the fluctuations. The 
backreaction of the growing modes on the background 
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fields is expected to shut the growth down at some point, 
but exactly when? Backreaction is also the only hope 
to make models which exhibit parametric amplification 
of super-Hubble cosmological perturbations compatible 
with the Cosmic Background Explorer (COBE) normal- 
ization [ljj. 

In this paper we will investigate the effects of backreac- 
tion on the growth of matter and metric fluctuations us- 
ing the Bassett and Viniegra model (0) as our toy model. 
We will study the growth of scalar field and scalar metric 
perturbations, including the effect of backreaction in the 
Hartree approximation. We carefully treat the evolution 
during inflation, which can be very important for super- 
Hubble scales. We will compare the large-scale normal- 
ization predicted for this model with the COBE value, 
and find that, although backreaction is crucial in limit- 
ing the growth of the fluctuations, the final amplitude 
is larger than allowed by the COBE normalization (for 
supersymmctry- motivated coupling constant values), un- 
less the period of inflation is very long. Note that the final 
amplitude of fluctuations in our model is independent of 
the scalar field coupling constant (unlike what happens 
without parametric resonance effects), but that it may 
depend on the duration of the inflationary period. We 
also extend the model to study the effect of x-field self- 
coupling, which can be important in limiting the growth 
of fluctuations. 



II. MODEL AND LINEARIZED DYNAMICS 

A. Equations of motion and analytical theory 

Our model is the two-real-scalar-field, gravitationally 
minimally coupled model specified by the Lagrangian 
density 
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The field <j) drives inflation, while \ is significant only af- 
ter parametric resonance begins, so the inflationary dy- 
namics is essentially that of (X/4)F chaotic inflation. 
Note that the behaviour of this system is expected to be 
robust under the addition of a small mass term ni^tfi 2 

with rrirj, <C %/Ampi and for the ratio of coupling con- 
stants satisfying g/V\ < \f\m-p\/m ( f, |lS| |. In particular, 
this will be the case for supersymmetric models, which 
motivate the choice g 2 — 2 A [Bof. In addition, we will 
show that large values of g/VAare in fact inconsistent 
with the significant amplification of super-Hubble modes. 
On the other hand, for g/V~X > -v/Ampi/m^, the theory 
of "stochastic resonance" for a massive inflaton may need 
to be applied 0]. 



It is traditional to separate the inflaton field into a ho- 
mogeneous, "classical" background 4>(t) and a perturba- 
tion 8(j)(x,t), which begins as sub-Hubble quantum vac- 
uum fluctuations early in inflation. It will be useful to 
make a similar separation for the field \, although we 
have no reason to expect a non-zero initial homogeneous 
X component in this model. The equations of motion for 
the homogeneous parts of the inflaton and \ fields are 
the Klein- Gordon equations, 



4> + 3ff0 + \<P 3 + g 2 X 2 <P = 0, 



X + ZH X + g 2 Fx = 0, 



(3) 
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with Hubble parameter H = a/ a, where a is the scale 
factor. To complete the background dynamics we must 
specify the evolution of the background spacetime metric. 
We assume a spatially flat Friedmann-Robertson- Walker 
universe, for which the 0-0 Einstein equation gives the 
Friedmann equation 
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For 4> ~ mpi, the universe undergoes slow-roll infla- 
tion, with H approximately constant and the scale factor 
a increasing approximately exponentially with time. As 
slow-roll ends, the "damping" term 3Hcb becomes less im- 
portant in Eq. (|jj|) and the field begins to oscillate about 
= 0. This marks the start of the preheating period. 
Averaged over several oscillations, the equation of state 
(in the absence of backreaction) is very nearly that of a 
radiation-dominated universe [H, and the amplitude of 
the inflaton's oscillations decays as a -1 . This is a conse- 
quence of the (near) conformal invariance of this mass- 
less model, which considerably simplifies the treatment 
of parametric resonance as compared with the massive 
case 

In writing the linearized equations of motion for per- 
turbations about the background, we will use the longi- 
tudinal gauge. For this model the metric can be written 



ds 2 = (1 - 2<S>)dt 2 - a 2 (t)(l + 2<f>)dxidx\ (6) 

so scalar metric perturbations are described by the single 
variable <I>. The momentum-space first-order perturbed 
Einstein and Klein-Gordon equations are 
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where <j>\ = <f>, <j>2 = X? = dV/d<pi, and comoving 
momentum k subscripts have been suppressed for clarity. 
Equations (R) and can be combined to give 
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which fixes $ once the matter fields are known. 

An important quantity in the study of the linear evo- 
lution of metric perturbations is the Bardeen parameter 
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which for fc/a <C -ff and single field models satisfies the 
"conservation law" M 
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where w = Pj p is the equation of state (p and P denoting 
energy density and pressure, respectively). When <&k can 
be neglected, Eqs. ([ll]) and ([l2]) can be combined to give 
the familiar result that the change in on super-Hubble 
scales over some interval of time is determined solely by 
the change in equation of state. 

In the absence of metric perturbations, the linearized 
dynamics in the model described above is known to ex- 
hibit parametric resonance during preheating Jl8| . To see 
this, it helps to take advantage of the near conformal in- 
variance of the model and rewrite the equations in terms 
of conformally scaled fields (pi = a(j>i and a dimension- 
less conformal time x = -v/A^o 7 ^ where rj = j dt/a and 
4>o is the amplitude of inflaton oscillations at the start 
of preheating. Then, with the x background and metric 
perturbations set to zero, the inflaton background and 
perturbed field equations (||) and (|^) become 
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where 



/c 2 /(A</>q) is a dimensionless comoving mo- 



mentum and primes denote derivatives with respect to 
the scaled conformal time x. Here we have ignored 
terms proportional to a" /a since preheating is a nearly- 
radiation-dominated phase in this model |3|. The con- 
formal field (p then undergoes constant amplitude elliptic 



cosine oscillations, while the perturbation equations are 
Lame equations [[22|, which are known to exhibit reso- 
nance within certain bands in parameter space 111]. In 
particular, the Sxk equation exhibits strong resonance for 
a band that includes k = for the supersymmetric point 
g 2 /X = 2, and weak "narrow" resonance in small-scale 
bands. On the other hand, 5<f>k exhibits narrow resonance 
for a sub-Hubblc momentum range, independent of the 
coupling constants. For resonant modes, the growth is a 
modulated exponential, Sxk °c e tikX , with Floquet index 
/ifc. If we allow g 2 /X to vary, we find a sequence of Sxk 
resonance bands for k = 0, centred at g 2 /A = 2n 2 with 
width 2n, for positive integral n jl8). The Floquet index 
reaches a maximum value of ^t max — 0.238 at the centre 
of each k = band. 



B. Numerical results 

For our numerical calculations, we were primarily in- 
terested in the behaviour of cosmological-scale matter 
and metric modes. Thus we evolved a scale which left 
the Hubble radius (at time to) at about N — 50 e-folds 
before the end of inflation. For (A/4)0 4 models, the num- 
ber of e-folds during slow-roll inflation after initial time 

t is pa 
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thus we used the homogeneous inflaton initial value of 
(f>(to) = 4mpi. We began the calculations with the modes 
still somewhat inside the Hubble radius, so the initial 
conditions for the matter field fluctuations were simply 
given by the conformal vacuum state 
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*&fr(*o) = -ivi(to)S<f>i(t ), 
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with lu 2 (t) = (k/a) 2 + 3A0 2 + g 2 x 2 and uj\ (i) = (fc/a) 2 + 

g 2 4> 2 . Physically, the a -3 / 2 dependence arises because 
particle number densities n k oc |<50ifc| 2 must decay like 
a~ 3 in the massive, adiabatic regime. The initial metric 
perturbations were then determined by Eq. (|l0|). 

To illustrate the dynamics in the absence of backreac- 
tion, we numerically integrated the coupled set of back- 
ground equations (||) and ^ and perturbation equations 
(|8|) and (|9j) using the initial conditions described above, 
and for g 2 /A — 2, A = 10~ 14 , and a zero x background. 
We used the constraint Eq. ( |l0| ) as well as the conserva- 
tion equation ( 
tions. In Fig 



to check the accuracy of the calcula- 
Ij we display the evolution of our cosmo- 
logical modes, together with the Bardeen parameter 
during inflation and preheating. For each of the pertur- 
bations Xk = Sxk, S4>k, &k, and ( k we plot the power 
spectrum f24] 
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FIG. 1. Numerical simulation of linear cosmological-scale 
perturbations in the two field model described in the text, 
in the absence of backreaction. Plotted are the logarithms 
of the power spectra Px(k) = (k 3 /2Tr 2 )\Xk\ 2 , for Xk — S\k, 
8<f>k, and ( k , and using mpi = 1, g 2 /X = 2, A = 10~ 14 , 
and zero x background. The main figure shows the evolution 
during inflation, while the inset details the behaviour during 
preheating, using the rescaled conformal time x. This partic- 
ular comoving scale leaves the Hubble radius approximately 
5 e-folds after the start of the simulation. 
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rather than the mode amplitudes, to facilitate compari- 
son with the COBE measured normalization which gives 

v* ~ i(r 10 £§]. 

Figure |l| shows how the modes begin early in infla- 
tion as sub-Hubble oscillations, and become "frozen in" 
after they exit the Hubble radius. Note that the 5\k 
fluctuation experiences some damping late in inflation, 
when its effective mass squared g 2 <fi 2 becomes somewhat 
greater than H 2 , which decreases like 4 in the slow-roll 
approximation. The inflaton perturbation 8<t>k, however, 
stays roughly constant during inflation even though its 
effective mass is comparable to that of the Sxk mode. 
This is because of the coupling between 5(f>k and 
in the linearized perturbation equations. We can also 
observe a growth of between the time the mode 
exits the Hubble radius and the beginning of preheat- 
ing, by a factor of approximately 20, in good agreement 
with the growth predicted from the "conservation law" 
Eq. Jl2]). Also note that after this small growth stage 
the cosmological-scale metric power spectrum ends up 
close to the 10~ 10 (~ e~ 23 ) level, as the standard the- 
ory predicts for A ~ 10~ 14 in the absence of parametric 
resonance M . During preheating we observe exponential 



60 



40 



20 



- 



-20 - 



-40 



-60 




200 



FIG. 2. Simulation of linear cosmological-scale perturba- 
tions in the absence of backreaction as in Fig. [j], but with 
a non-zero initial background x(^o) = 10 _10 mpi. Plotted are 
the logarithms of the power spectra Vx{k) = {k 3 /2n 2 )\X k \ 2 , 
for X k = &Xk, 8<t>k, and <E> fc , and using mpi = 1, g 2 /\ = 2, 
and A = 1CP 14 . Coupling through the x background causes 
6<f>k and <E>fc to grow initially at twice the rate of Sxk- 



growth of Sxk while the super-Hubble S(f>k mode does not 
grow, as expected from the analytical theory. and 
remain constant, since according to Eq. (|8|) the metric 
perturbations couple only to 8<j)k in the absence of a \ 
background, at linear level Jl6| . 

To observe the effect of including a non-zero homoge- 
neous x background, we repeated the above calculation 
with an initial value of x(io) = 10 _10 mpi (this value il- 
lustrates well the various stages of evolution). Figure 
§ indicates that Sxk grows as before, but 5<fik and 
now grow initially with twice the Floquet index of Sxk- 
This is the result of the driving term 2g 2 (j>x^Xk in the 
equation of motion for 5<f>k , Eq. (|9|) , which contains two 
factors growing like e Mmax:r (clearly the evolution of the 
background x w iU be essentially the same as that of 5x k 
for k <C aH). Once the background x field becomes com- 
parable to the inflaton background, all the perturbations 
synchronize and grow at the same rate. We will discuss 
the significance of the homogeneous x field in relation to 
the nonlinear evolution of the fields in the next section. 



III. BACKREACTION 

A. Equations of motion 

The linearized equations in the previous section de- 
scribe unbounded growth of perturbations during res- 
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onance. In reality this growth must of course stop at 
some point, namely when the perturbed field values are 
on the order of the background values. A full nonlin- 
ear simulation will include this effect automatically, but 
approximation methods can alleviate the computational 
costs significantly. A common approach to approximate 
this backreaction on the background and perturbation 
evolution is to include Hartree terms in the equations 
of motion 25 |]. This entails making the replacements 
& -»• & + 8$ -> (6$), and <50? -> 3(Stf)5fa. In 
this approximation, the background equations (||) - (|^) 
become 
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^W) + ^x 2 < 



(20) 



+ 3fT0 + + 3A(<50 2 )0 + .g 2 (<Sx 2 )0 = 0, (21) 



X + 3Hx + V x +g 2 (5<b 2 )x = 



(22) 



Similarly, the momentum-space linearized field perturba- 
tion equations @ become 

+ V ,M S< t>ik = 40$ fc - 2y $ fe , (23) 

3 



<5 Xfe + 3H8 Xk +{-2+ <TW> ) Sxk 



(24) 



In this approach, the field fluctuations are calculated self- 
consistently from the relations 
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In practice, the resonance band will provide a natural 
ultraviolet cutoff. 

The Hartree terms approximate the full nonlinear dy- 
namics of the fields. To illustrate what this approxima- 
tion entails, we may consider the exact dynamics of the 
fields, treating the Klein-Gordon equation as a classical 
field equation. This should be a good approximation soon 
after the beginning of the resonance stage, since occupa- 
tion numbers will grow exponentially with time p|,|l8|]. 
As an example, consider the exact evolution equation for 



64> in position space, obtained by perturbing the Klein- 
Gordon equation, setting the background X to zero, and 
ignoring metric perturbations, 



3H5(f> 
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V 2 <5(/> + 3X(j) 2 



3\4>6<p 2 + XScj) 3 



+ g 2 5 X 2 + 9 2 8 X 2 5<P = 0. 



(26) 



The terms in this equation describing the interaction be- 
tween the <j> and x fields become in momentum space 
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Thus the Hartree term g 2 {8x 2 )$4 l k in Eq. ( p3[ ) corre- 
sponds to the second term in expression (p7|), restricted 
to k" = — k'. Physically, this means that only scattering 
events which do not change the 8(j) k momentum are in- 
cluded in the Hartree approximation, and "rescattering" 
events are ignored. 

It is important to notice that the first term in (|27 
which scatters particles from the homogeneous inflaton 
background into mode 8<f> k , could be larger than the 
Hartree term since initially \<f>\ > \8<fi\, unless the first 
term vanishes upon averaging (integrating) over the en- 
tire phase space of contributing terms (which is what is 
assumed in the Hartree approximation). If it does not 
vanish, the first term in ( p7j ) will act as a driving term 
for the 8(f) modes in (p6|). Since some 8x modes experi- 
ence parametric amplification with Floquet exponent /x, 
this term will lead to an important second-order effect, 
namely the growth of 8<f> as e 2 ^ x . This effect is left out in 
the Hartree approximation. Because the metric pertur- 
bations are coupled to 5<j) through Eq. (0), we also expect 
that, with the homogeneous X se t to zero, the Hartree 
approximation will miss the corresponding growth of 
Note, however, that by including the X background term 
2g 2 x4>3x m (p6|); and setting % 2 ~ (<5x 2 ), we can approx- 
imate the effect of the important first term in fl27| ) , as we 
saw in Fig. g. 

Just as with the scalar fields, metric fluctuations may 
grow rapidly in our model. We can account for the 
backreaction of metric perturbations through the effec- 
tive energy-momentum tensor formalism of Abramo et al. 
p6| . This involves expanding the Einstein equations to 
second order in the perturbations and taking the spa- 
tial average to obtain effective background equations. In 
our case, the metric and inflaton equations (^0|) and ( |2~l| ) 
become, with background \ set to zero, 



H 2 = 
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■ 3Hj>)(l + 4($ 2 )) + X(t) 3 + 3\(5<j) 2 )4> 

■ 9 2 (<>X 2 )(t> ~ 2($<50) - 4{$6<t>) 

6H{$6</>) + 40($$) - -^(®V 2 6(j)) = 0. (29) 



B. Analytical estimates 

1. Evolution of perturbations during inflation 

Perturbations will grow during parametric resonance 
until backreaction becomes important. We can analyt- 
ically estimate the amount of growth by estimating the 
time at which the Hartree term g 2 (Sx 2 ) is of the order of 
the background A</> 2 (cf Eq. (|23|)). Note that in the ab- 
sence of metric fluctuations, such an estimate should be 
accurate, at least for q 2 /A ~ 1, as nonlinear lattice simu- 
lations indicate |27|2^. In addition, we expect the matter 
sector to dominate the dynamics. In order to estimate 
the variance (Sx 2 ), we wm need to calculate the evolu- 
tion of Sxk modes, starting from the adiabatic vacuum 
inside the Hubble radius, continuing through inflation, 
and finally through preheating. The evolution during in- 
flation is quite complicated, and will have a crucial effect 
on the final variances, so we will describe the inflationary 
stage in some detail. We consider general values of g 2 / A, 
rather than just the supersymmetric point. 

We will only need to consider the contribution to (Sx 2 ) 
from modes which are super-Hubble at the start of pre- 
heating. To see this, first note that for .g 2 /A ~ 2, the 
small-scale boundary of the strongest (and largest-scale) 
resonance band is at fc max /i — vA0o/2, where (f>o{t) is 
the amplitude of inflaton oscillations during preheating 
p8| . Next, we can use the Friedmann equation (||) to 
write the Hubble parameter in terms of 4>o, giving 



H 2 = 

3mpj 



(30) 



(Note that this equation also applies approximately dur- 
ing slow-roll.) Using the value (f>o = 0.2mpi, we calculate 
the ratio ai//fc max ~ 0.6 at the start of preheating. Thus 
the Hubble radius corresponds closely to the smallest res- 
onant scale. This result is not very sensitive to g 2 /X as 
long as we are near the centre of a band, i.e. g 2 /X = 2n 2 , 
since fc max increases only slowly with g 2 /X in this case 
jf8|. Also, we can ignore the resonance bands at higher 
k values, since they correspond to narrow resonance. 

To estimate the evolution of 8\k on super-Hubble 
scales during inflation, we can ignore terms containing 
the background x as wen as the spatial gradient term in 
Eq. (|^), resulting in a damped harmonic oscillator equa- 
tion with time-dependent coefficients, 



5x k + 3H5 Xk + g 2 2 Sxk = O 



(31) 



During slow-roll, we can use the adiabatic approximation 
to find solutions to this equation, since \H\ <C H 2 . Thus 



for g 2 <\> 2 > (3H/2) 2 we have underdamped oscillations 
with damping envelope 



Sxk oc exp 



(3H/2)dt 



-3/2 



(32) 



For g 2 (j) 2 < (3H/2) 2 , we have the overdamped case with 
two decaying modes. Ignoring the more rapidly decaying 
mode, we obtain 



Sxk oc exp 



J (3H/2 - V9ff 2 /4-.g 2 2 ) dt 



(33) 



In this case, the fluctuations are very slowly decaying in 
the massless limit g 2 <p 2 <C (3H/2) 2 , while they approach 
the a~ 3 / 2 decay as g 2 <\> 2 -> (3H/2) 2 . 

During slow-roll we have H 2 oc (j> 4 (see Eq. (|30|)), so 
that H 2 decreases more rapidly than g 2 (f> 2 , and there is 
a transition between the over- and underdamped stages. 
The two types of behaviour are separated by the critically 
damped case, g 2 <j) 2 = (3H/2) 2 . Using Eqs. © and (||), 
we can write this critical damping condition in terms of 
the number of e-folds after critical damping, N cr n, as 



N clit = In 
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^crit 



3 A \ k cr i t 



(34) 



where subscript "f" refers to the end of inflation and 
"crit" to the time of critical damping. Wavevectors fc cr it 
and kf leave the Hubble radius at i cr it and tf, respec- 
tively. We see that as <? 2 /A increases, cosmological scales 
are damped like a" 3 / 2 during a greater and greater part 
of inflation. We thus expect that for large enough g 2 /X, 
the backreaction of the smaller-scale modes will termi- 
nate parametric resonance when cosmological-scale Sxk 
modes are still greatly suppressed. In other words, there 
will be a maximum value of g 2 /X for which there is sig- 
nificant amplification of super-Hubble Sxk perturbations, 
as anticipated in |rij| . 

We first consider the evolution of the modes which 
leave the Hubble radius after t cr it, *.e. k > k CT i t (but 
which are still super-Hubble at the end of inflation, 
k < kf). These modes are effectively massive during infla- 
tion, and hence we can simply use the adiabatic vacuum 
state, Eq. (17), which for k <g; aH gives 
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Note that if we define the spectral index n through 
V x {k) oc [Q, then for this part of the spectrum 

we have n = 4, an extreme blue tilt. 

Next, we will calculate the evolution of modes which 
leave the Hubble radius before t cr it, i.e. modes with k < 
fc cr ; t . In this case, the modes are approximately massless 
when they exit the Hubble radius (g 2 (j) 2 < (3H/2) 2 for 
t < tcrit), so we can use the standard result for a massless 
inflaton 



G 



\5xk(t k f 



H 2 (t k ) 
2k 3 : 



(36) 



where t/. is the time that mode S\k exits the Hubble 
radius. We now must use Eq. ([33|) to evolve the modes 
during the overdamped period, % < t < t alt . Writing 
dt = d<f>/4>, and using the slow-roll approximation cf> ~ 
—V,4>/3H, we can perform the integral to obtain 



H 2 (t k ) 
2k 3 ' 



-3F(N k ) 



(37) 



where N k is the number of e-folds after time t k and 
F(N k ) = N k - N clit - ^/N~ k y/N k ~ N crit 



+ N crit In 



{Nk + VN k - N c , it 
VN^Tt 



(38) 



Next we can readily propagate the modes through the 
underdamped period, f cr it < t < if, using Eqs. ( |32| ) and 
(pi), giving 



\SXk(tt)\ 



2 _ H g-3F 



(N k )-2g 2 /\ 



2k 3 



(39) 



Since the damping term F(N k ) is positive, we see as ex- 
pected that large-scale modes are strongly damped for 
large g 2 /A. 

Finally, we can approximate the conformal time de- 
pendence of all super-Hubble modes during parametric 
resonance as 



(40) 



if we are near the centre of a resonance band. This is 
valid since, in this case, the Floquet index [i k varies only 
slightly for scales larger than a few times the Hubble 
radius (i.e. the smallest resonant scale) fl8| . 



2. Variances and total resonant growth 

Now we can proceed to calculate the field variance, 
(Sx 2 )- We will use Eq. (Bp), restricting the integral to 
the resonantly growing modes. We begin with the case 
g 2 /X = 2. Equation ( |34| ) tells us that in this case 
N CI it ~ 4/3, so that essentially all of the evolution during 
inflation is in the overdamped regime, and we only need 
to consider modes with k < fc cr it- The variance integral 
will be dominated by modes with N k 3> N CT i t , so we may 
approximate the damping term in Eq. ( 38 ) as 



-3F(N k 



Nr. 



4Mt 



s 2 M 



(41) 



For the current case, g 2 /A = 2, we can now combine the 
expression © with Eqs.©, ©, ©, and © to 
obtain for the power spectrum on resonant scales at the 
end of preheating 



^(Me) = ^e 2 " • (42) 

Here t c is the time that the resonance shuts down, and 
x c is the corresponding scaled conformal time. As we 
will see, the important thing about this result is that the 
power spectrum is essentially Harrison-Zel'dovich (inde- 
pendent of fc), with spectral index n = 1. 

We can next rewrite the variance integral, Eq. (p5|), in 
terms of the power spectrum as 



(Sx 2 (t e )} 



N„ 



dN k V x (k,t c ) = N V x (t e ), 



(43) 



where N ~ 50 is the total number of e-folds during 
inflation. Finally, the criterion g 2 (Sx 2 (t c )) ~ \<f> 2 (t c ) 
gives, using the value 4>{t e ) 



10- 2 m P i, 



10" 6 m 2 



pi 



(44) 



for the $Xk power spectrum on cosmological scales at the 
end of preheating. Note that this result used only the 
fc-independence of the power spectrum (which is a result 
of the special choice g 2 /\ = 2), and the values of iVo 
and 4>{t e ). In particular, the result is independent of A, 
unless, contrary to our implicit assumption, A is so large 
that g 2 (Sx 2 ) > A0 2 already at the start of preheating. In 
this case, Eq. (E3) will be an underestimate. 



According to the results from Section II B, we expect 
synchronization of the other fields to 8xk , so that in par- 
ticular we expect P$ ~ V x /m,py Therefore we conclude 
that, for g 2 /X = 2, the metric perturbation amplitude 
will indeed be considerably larger than the COBE mea- 
sured value, even including the effect of backreaction. 

Next we will repeat the preceding analysis for the sec- 
ond super-Hubble resonance band, at ,g 2 /A = 8. In this 
case we have iV cr it ~ 5, so we must consider modes that 
exit the Hubble radius both before and after t cr j t . For 
the large-scale modes, k < fc cr it> it will be sufficient to 
place an upper limit on the variance. Using Eq. ([39]) , but 
ignoring the damping factor e~ 3F , we obtain 



V (k tA < g2 ( fCTit ) r -2fl 2 /A 

xl ' f) (2tt) 2 



2 x 10" 8 Am| 1 (45) 



on scales k < fc cr j t at the end of inflation. Thus, using 
Eq. d43|), the contribution to the variance from modes 
with k < fccrit satisfies the (probably very conservative) 
bound 



(Sx 2 (ti))k<k CIit < 9 x 10- 7 Ato 



(46) 



Next we can use Eq. (|35|) to calculate the contribution to 
(<5% 2 ) from smaller-scale modes with fc crrt < k < kf, 



(^X 2 ( i f))fc cri t<fc<fc t 



167T 3 'a 3 gw Jkci . t 




d A k 



3 x 10~ 6 Am 2 



(47) 

(48) 
(49) 
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Here we have used fcj? rit <C kf (which follows from Eq. j3^ ) 
for g 2 /X = 8), the relation fc f /a f = H(t{), Eq. (|30|), and 
the value 4>f = 0.2mp\. This value of the small-scale vari- 
ance exceeds our upper limit on the large-scale variance 
in Eq. ( [l6| ) , so we can ignore the contribution from the 
large-scale modes, {Sx 2 {tf))k<k crit ■ 

Now we can again apply the condition g 2 (5x 2 {t c )) ~ 
X(j} 2 (t c ), which in this case gives 



4A" 



(50) 



Finally, we can use Eq. ( p9|) without approximation to 
calculate the cosmological-scale power spectrum at the 
end of preheating, for the case g 2 /X = 8, 



V x (t e ) 



H 2 (t Q ) 
(2n) 2 



■ exp 



'pi- 



-3F(N ) - 2y + 2 Mn 



(51) 
(52) 



In this case the growth stops before the cosmological per- 
turbations exceed the COBE value, and thus paramet- 
ric resonance does not change the standard predictions 
j?J for the size of the fluctuations. Therefore, since the 
damping of super-Hubble 5\k modes increases as g 2 /X 
increases, the standard predictions are not modified for 
all resonance bands beyond the first, i.e. for g 2 /X > 8. 

To close this section, consider the behaviour of the 
large-scale variance if we suppose that inflation started 
much earlier than the time that cosmological scales left 
the Hubble radius, i.e. let Nq — > oo. In this limit we see 
from the approximation Eq. ([fl|), the expression H(tk) oc 
Nk, and Eq. (B9|) that the variance becomes 



(Sx Z )k<k crlt oc / dN k N< 



N<> 



(53) 



which, for g 2 / X = 2, is divergent as Nq — > oo. Such 
divergences are well-known in inflationary models |?0| . 
Here the divergence suggests that for large enough A^o, 
the growth of cosmological-scale modes will stop before 
they exceed the COBE amplitude, due to the large con- 
tribution to the variance from super-cosmological scales. 
Indeed, for No ~ 10 6 (a value not out of the question in 
chaotic inflation Q) Eq. @ gives V x (t c ) ~ 10 -10 ro£i 
for g 2 /X = 2. For g 2 /X > 3 the variance converges, al- 
though in this case cosmological-scale modes are already 
supressed by the mechanism described above. 



C. Numerical Results 



we begin each mode's evolution inside the Hubble radius 
during inflation, using the initial vacuum state, Eqs. Jl7| ) 
and (|lS|). Each mode is incorporated into the calcula- 
tion shortly before it leaves the Hubble radius, so that 
the spatial gradient terms are never too large. The vari- 
ances are calculated by performing the discretized inte- 
grals, Eqs. (p5|), only over the resonance band; thus they 
are convergent. Note that the variances are calculated 
simultaneously with the field backgrounds and perturba- 
tions. 

In Fig. ^| we present the evolution of the 5\k, o~4>k> 
and $fc power spectra on the same cosmological scale as 
was studied in Section II. All parameters are the same 
as for Fig. ||, except here we use for the initial back- 
ground value x(to) — 10 _6 mpi, which means that during 
preheating \ 2 ~ {&X 2 )- The evolution is initially simi- 
lar to that of Fig. g, only here the growth saturates at 
in good agreement with our predic- 



3 x 10 



tion based on Eq. (|4J). Also, as expected, the other fields 
closely follow V x . Whereas in the linear calculations 
the Einstein constraint equation ( |l0| ) was satisfied to ex- 
tremely good accuracy, with the inclusion of backreac- 
tion V<s> saturates at a factor of roughly 10 3 higher using 
Eq. ( |l0| ) than the illustrated result, which used Eq. (^). 
Note that a similar result was found in ]l3| ]. We sus- 
pect that this is a fundamental problem related to our 
attempt to capture some of the nonlinear dynamics with 
the Hartree approximation. Regardless of which value is 
used, the cosmological metric perturbations considerably 
exceed the COBE normalisation. In addition, we find no 
significant difference in the results when backreaction of 
metric perturbations is included using Eqs. ( p8| ) and j2S|), 
as expected if the matter fields dominate the backreac- 
tion. Thus all of our presented results exclude the metric 
backreaction terms. 

As discussed above, larger values of g 2 /X result in in- 
creased damping of Sxk on large scales during inflation, 
and at large enough g 2 /X we expect insignificant ampli- 
fication of super-Hubble modes. This is illustrated in 
Fig. U Here we examine the second resonance band at 
g 2 /X = 8, but use otherwise identical parameters to Fig. 
||. Resonance stops at V x ~ 10 _14 TOp[, consistent with 
our analytical estimate from Eq. (|52"|), and not exceeding 
the standard predictions for A ~ 10~ 14 0. Note that the 
small rise in at late times should not be trusted, as 
our Hartree approximation scheme will not capture the 
full nonlinear behaviour. For resonance bands at even 
higher g 2 /A, we find extremely suppressed cosmological 
Sxk amplitudes, in quantitative agreement with the cal- 
culations of the previous section. 



It is straightforward to check our analytical estimates 
from the previous section by numerically integrating the 
coupled set of Hartree approximation evolution equations 
( pp| ) - (^5|) and metric perturbation equation (||). We 
now must evolve a set of modes that fill the relevant 
resonance band. For example, for g 2 /X = 2, the first res- 
onance band extends from k = to n = 0.5 M]. Again 



8 



IV. SELF-INTERACTING X MODELS 
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FIG. 3. Numerical simulation of cosmological-scale pertur- 
bations with Hartree backreaction terms included and with 
a non-zero initial background, x{to) = 10 _6 mpi. Plotted are 
the logarithms of the power spectra Vx{k) = (k 3 /2TT 2 )\X k \ 2 , 
for X k = Sxk, Stj) k , and $ fe , using topi = 1, g 2 /\ = 2, and 
A = 10~ 14 . Backreaction terminates the growth of each field 
perturbation. 
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FIG 4. Same as Fig. 



in the second resonance band, 
X(t ) = 10~ 6 mpi and A = KT 1 



except for parameters lying 
. e. g 2 /\ — 8, and with 
Here backreaction of the 



small-scale modes terminates the growth before the COBE 
value is exceeded. 



A. Positive coupling 

We now consider the addition of a quartic self- 
interaction term for the x held, so that our potential 
becomes 



A 



A 



V, 



(54) 



with g 2 > 0. The significance of such a term for para- 
metric resonance was studied in lattice simulations j2j| 
and analytically |3l| l, but in the absence of metric per- 
turbations. Bassett and Viniegra [JL7| included metric 
perturbations, but ignored backreaction. Essentially, for 
A x > A we expect the x self-interaction to limit the 



growth of perturbations as compared with the A x 







case studied above, due to the presence of the "potential 
wall" (A x /4) x 4 . 

More precisely, the linearized equation of motion for 
the x field perturbation becomes, with x self-interaction 
but ignoring metric perturbations, 



5 Xk + 3H6 Xk +[- + 3X xX z + g z P ) S X k 



25 



= 0. 



(55) 



Thus for small enough initial x background, the initial 
behaviour of the modes will be essentially unchanged 
from the A x = case. However, when the x background 
grows to the point that x 2 /4> 2 ~ ff 2 /A y , the analytical 
parametric resonance theory of Section II A no longer 



applies, and we may expect the perturbations to stop 
growing. Since, as discussed above, for the significant 
production of supcr-Hubblc modes we require g 2 ~ A, 
we expect that x self-interaction will shut down the reso- 
nance when x 2 /^ 2 ~ A/A x , as long as X x > A. If A x < A, 
then the % 4 interaction term will not lead to a shutdown 
of the resonance since (based on our numerical simula- 
tions) the homogeneous x field never substantially ex- 
ceeds the value of the inflaton background. 

We have confirmed this expectation numerically, and 
we give an example of our results in Fig. ^. Here we have 
included Hartree backreaction and metric perturbations, 
and used coupling constant values A = 10~ 14 , .g 2 /A = 2, 
and A x = 10~ 10 , and initial backgrounds </>(to) = 4mpi 
and x(^o) = 10~ 6 mpi. We indeed observe the termina- 
tion of the super-Hubble modes' growth at approximately 
the time when x 2 /0 2 = A/A x . 

Note that, in the absence of backreaction, Bassett and 
Viniegra observed a continued slow growth of super- 
Hubble perturbations after the initial termination of the 
resonance when x 2 l<t> 2 ~ A/A x p7| . We confirmed this 
result; however, we note that when we include the back- 
reaction term 3A x (<5x 2 ) in the evolution equations, we ex- 
pect backreaction to become important also at the time 
that x 2 /^ 2 ~ A/A x , with our choice x 2 — i^X 2 )- Hence, 
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FIG. 5. Same as for Fig. |3| except with x field self-coupling 
A x = 10 -10 . The x self-coupling causes the early termination 
of growth. 



as seen in Fig. pL the slow growth is completely sup- 
pressed. 



B. Negative coupling 

The presence of x self-coupling means that we no 
longer require g 2 > for global stability. In fact, for 
the case g 2 < 0, the potential will be bounded from be- 
low for AA X / g 4 > 1 . This negative coupling case was 
studied in the absence of metric perturbations using lat- 
tice simulations in f32|| , and without backreaction in [fl6[ . 
The behaviour of the fields is qualitatively different in 
the negative and positive coupling cases. For g 2 < 0, 
potential minima exist with non-zero homogeneous part 
of the x field. Thus, assuming the fields fall into these 
minima, the problem of choice of x background discussed 
in previous sections for the positive coupling case is alle- 
viated. 

For initial homogeneous \ fields large enough (x(to) ^ 
mpi), we find numerically that the fields fall into the 
potential minimum by the end of inflation, and the 
two fields subsequently evolve in step during preheating. 
This effectively reduces the system to a single-field sys- 
tem, and hence no resonance is possible on super-Hubble 
scales. 

To see this explicitly consider the case A x = A, for 
which the symmetry of the potential requires the po- 
tential minima to lie along x 2 = 4> 2 ■ If we choose the 
same initial signs for x an d 4>i then during preheating 
the backgrounds lie along the attractor <fi = x- Similarly, 



since the behaviour of super-Hubble modes is essentially 
the same as that of the backgrounds, we have Sx = §<t> 
during preheating. Then the perturbation equation ( |55| ) 
becomes 



ri.2 



Sx k + 3ff(5x fc 



+ 3(A + g 2 



6 X k = 0. (56) 



Thus the effective mass of the Sx oscillations is precisely 
three times the effective mass of the background inflaton 
oscillations (cf Eq. (||)), so that just as with the case 
of the inflaton perturbations in Eq. (|lj) , there will be 
no resonance on super-Hubble scales for all allowed val- 
ues of g 2 . We have confirmed this numerically; indeed 
more generally, as long as initially x(^o) ~ m P\ but for 
any A x > A, the two fields will be proportional during 
preheating and no super-Hubble resonance will result. 

This result assumes that during preheating only the 
"field" Sx + S(f> is excited. If orthogonal field excitations 
Sx — S(f> are present, they can grow resonantly. The effec- 
tive squared mass of Sx — S(j) excitations is (3A — g 2 )(j) 2 , 
so that according to th e analytical parametric resonance 
theory of Section II A, super-Hubble resonance will oc- 
cur near 3 A — g 2 = 2n 2 (X + g 2 ), for integral n (we require 
n > 2 for negative g 2 ). That is, super-Hubble Sx — S<f) 
modes will grow for g 2 ~ A(3 — 2n 2 )/(l + 2n 2 ). However, 
numerically we observe only extremely small components 
Sx — S<f> by the end of inflation, so their growth is sub- 
stantially delayed. 

On the other hand, for small initial homogeneous part 
x(*o) ^ mph wc nn d that the potential minima are 
not reached by the end of inflation, and the two fields 
evolve in a very complicated manner during preheat- 
ing. The analytical theory of parametric resonance can- 
not be applied, but numerically we do find roughly ex- 
ponential growth of super-Hubble modes in this case, 
as found in [pUf . The growth rate increases as g 2 de- 
creases towards the value at which global instability sets 
in, g 2 ~~ 



We have illustrated this case in Fig. p|, using the 
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-0.5A, A, 



A. 



parameter values A = 10 
</>(io) = 4mpi, and x(*o) = 10~ 6 mpi. Here the growth 
rates and final power spectra values arc comparable to 
the A x = case of Fig. ||, though the Sx field is 
not damped during inflation for negative coupling. For 
A x > A, the growth is terminated early, just as in the 
positive coupling case. 



V. SUMMARY AND DISCUSSION 

In this paper we have studied backreaction effects on 
the growth of super-Hubble cosmological fluctuations in a 
specific class of two field models with a massless inflaton 
4> coupled to a scalar field x- Otrr study was based on 
the Hartree approximation. 

For the non-self-coupled x field case, we found that 
backreaction has a crucial effect in determining the final 
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FIG. 6. Same as for Fig. ^, except for the negative coupling 
case, with parameters A = 10 -14 , g 2 = — 0.5A, and A x = A. 
Here the growth is comparable to the A x = case of Fig. ^| 
even though analytical parametric resonance theory cannot 
be applied here. 

amplitude of fluctuations after preheating. For values 
of the coupling constants satisfying g 2 /X = 2 (the ratio 
predicted in supersymmetric models), the predicted am- 
plitude of the super-Hubble metric perturbations at the 
end of preheating is too large to be consistent with the 
COBE normalization, thus apparently ruling out such 
models. One possible loophole is the backreaction con- 
tribution from super-cosmological scale fluctuations. For 
sufficiently long periods of inflation, the predicted am- 
plitude can be consistent with the COBE normalization. 
In addition, the final amplitude of the fluctuation spec- 
trum is independent of the coupling constant A. Note 
that the growth of inflaton fluctuations 5<j>k (and hence 
metric perturbations $fc) occurs in these models either 
through coupling to 5\k via a homogeneous background 
X field or through nonlinear evolution effects. 

The situation for g 2 / A 3> 1 is very similar to the pre- 
viously studied case of a massive inflaton in the broad 
resonance regime ]l2[-p^| . Cosmological-scale 8\k modes 
are significantly damped during inflation, and the end of 
resonant growth is determined by the growing small-scale 
modes. Already for the second resonance band (centred 
at g 2 /A — 8) cosmological metric perturbations are not 
amplified above the COBE normalization value. This 
implies that preheating does not alter the standard pre- 
dictions for the $k normalization in (A/4)<^ 4 inflation for 
the second and all higher resonance bands. The impor- 
tant difference between the model we have studied and 
the massive inflaton case is that, in the massive model, 
weak super-Hubble suppression at small g 2 is accompa- 



nied by weak resonant growth during preheating [T^ ] , so 
that no significant super-Hubble amplification is possible. 

The inclusion of x ne ld self-interaction alters the evo- 
lution in a predictable way: the resonant growth stops 
when x 2 /^ 2 ~ A/A x , as long as A x > A. This means that 
we are unable to rule out models (on the basis of a too 
large production of metric perturbations) with A x /A > 
10 4 . In the negative coupling case, there are two pos- 
sibilities. For large initial x backgrounds, x{to) ~ m ¥h 
the system becomes essentially single-field, and no reso- 
nance occurs (at least until late times). For small initial 
X, exponential growth occurs for large enough allowed 

\9 2 l 

The Hartree approximation provides a useful approach 
for the inclusion of the effects of backreaction. However, 
as mentioned in Section III, this approximation misses 
terms which could contribute to the evolution of fluctu- 
ations in an important way. We believe, nevertheless, 
that our results are sufficiently accurate to predict, for 
the models studied, whether or not the metric pertur- 
bation amplitude after preheating is consistent with the 
COBE measurement. Nonlinear effects, or rescattering, 
will primarily affect the detailed evolution of matter fields 
after backreaction is important. Still, it is of great inter- 
est to extend our analysis to a full nonlinear treatment, 
as was done in the absence of gravitational fluctuations 
in [ p3p^j2^ ], and including metric fluctuations in || for 
single field models. 
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